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Abstract 

Motzkin numbers are derived from a special case of Random Domino Automaton - recently pro- 
posed toy model of earthquakes [jj. An exact solution of the set of equations describing stationary 
state of Random Domino Automaton in "inverse-power" case is presented. A link with Motzkin 
numbers allows to present explicit form of asymptotic behaviour of the automaton. 

PACS numbers: 02.30.Lt (Sequences, series, and summability), 02. 10. Ox (Combinatorics; graph theory), 
02.50.Ey (Stochastic processes), 45.70.Ht (Avalanches), 

Keywords: stochastic cellular automata, Motzkin numbers, avalanches, solvable models, integer sequences, 
toy models of earthquakes 



* Electronic address: bialecki@igf.edu.pl 



It is known, how to obtain Catalan numbers out o: 



Domino Automaton 
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square lattice [2|. Here we derive Motzkin numbers js-G] from the recently proposed Random 



" the bond directed percolation on a 



8[, which may be regarded as a toy model of earthquakes. 
Random Domino Automaton comes from very simplified view of earthquakes. The space 
- one dimensional lattice - corresponds to boundaries of two tectonic plates moving with 
relative constant velocity. Due to irregularities of surfaces, relative motion can be locked at 
some places producing stress accumulation. Beyond some threshold of the stress, a relaxation 
took place. The size of relaxation depends on the nearby accumulated stress. Energy in the 
automaton is represented by balls added to the randomly chosen cell (each one is equally 
possible) with a constant rate - one ball in one time step. If the chosen cell is empty, it 
becomes occupied with probability v or the ball is scattered with probability (1 — z/) leaving 
the state of the automaton unchanged. If the chosen place is already occupied, there are 
also two possibilities: the ball is scattered with probability (1 — /x) or with probability [i 
the incoming ball triggers a relaxation - balls from the chosen cell and all adjacent occupied 
cells are removed. An example of relaxation of size five is presented in the diagram below. 



time = t 



time = t + 1 



The stationary state of the system may be described by the distribution of clusters. The 
number of clusters of the length for i = 1, 2, . . ., is denoted by nf, the number of empty 
clusters of length 1 is denoted by n^. Then the number of all clusters n and and the density 
p are 

i>l i>l 
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The following set of equations is derived from the stationarity conditions (see [l|]) 

"1 = ((1 - P)iV - 2n + n?) , (2) 



2 / n? \ 



1 

ni = T, X 



n / ^ — ' 

k=l 



for i > 3, where 



2n 

j2 = 

1 _^ 2Ei>iMi»"i 



From the above set the balance equation for the total number of clusters n and for the 
density p can be derived 

(1 - p)N ~2n = Y^ —riii, (5) 

^i'^- P) = j^C^f^inii^)- (6) 

i>l 

In the case which refers to equal probability of triggering an avalanche for each cluster, 
the parameters are fixed as follows: u = const and /i = |, where 6 = const. Equations ([6]) 
and dS]) give p = {0 + l)'^, and n = N9[{e + 1)(^ + 2)]-i, where ^ = ^ 6 (0, oo). Together 
with simple form of = 2n/(3 + 29) it allows to reduce the set of equations ©-(jl]) to the 
following recurrence: 

- ATg / 2 \ 

^1 - (^^ + i)(^ + 2)2^^+(2^^ + 3)y'' 
2 /2^ + l\ 

= TTV^ TTTT^ (8) 



9 + 2 \29 + 3 
2 [29 + 1 



+ 2) \29 + 3 
2 f 9 + l\^ 



fc=i 



for 2 > 2. 

We define new variables Cj for i = 0, 1, . . . , by 

/3 



^i+i, (10) 



where 



'2e + i\ , , 

« = 72^ > (11) 



+ 2) V2^ + 3 

B = + + (12) 

^ Ne{2e + i) ■ ^ ' 



Then, the equation (jH]) can be rewritten in the form 



m 



Cm+2 — Cra+l + / CkCm—kj (13) 

which is vahd for m > {m = i — 2) . Initial data Cq and Ci are easily obtained from equations 
(l7])-(j8]), when it is transformed according the rule of equation ( !T0|) . namely 



1 + 19 + 9' 

= = i + le + Aer ^''^ 

The above equation ( fT3l) has the form of Motzkin numbers recurrence SHa]. It is sim- 
ilar to the ubiquitous Catalan numbers recurrence (however the order of ffT3l) is not 1 

n 

but is 2) and can be solved using generating functions technique [3|. For the limit case, 
when ^ = 0, we start with cq = Ci = 1, and recurrence f|T3l) produces Motzkin numbers: 
1,1,2,4,9,21,51,127,323,835,2188,..., etc. Thus, it is explicitly shown, how to obtain 
them from the Random Domino Automaton. 

Below we present solution of the set ©-(IHl). The generating function 



m>0 



for the recurrence (fT3ll is equal to 



(1 - ^) - ^l-2z+{l- Aco)z^ + 4(co - ci)z^ 



2z' 

For Co = ci it reduces to 



(16) 



[1- z) - ^l-2z{l + kz) 
^2? 



G{z) = ^ ^ \. ' ^ (17) 



where k = (2co — |)- From expansion of the formula (ITB]) one can read explicitly the form 
of coefficients Cm being solutions of the recurrence ( |T3l) . Using 

(1 + ^)'' = EQ-"' (18) 

n>0 ^ ^ 



where 

a{a — l){a — 2) . . . [a — n + 1) 



for m > we have 



n / n\ 



(19) 



[z-]gf{z) = [."+2] l^-\^l-2z{l + kz)^ . (20) 



Then it follows 

/I 



1 - 2z{l + kz) = ^n>Q(^^ (-2)"2"(1 + kzY = 



n>2 

For any m > it gives 



Cm = [Z""^'] 



n>2 ^ ^ j=0 

Changing indices n + j = m + 2 for m > 1 we have 



r m+2 1 

2 J /o„ l\j 



(2c -1)^ /2{m-j) + l\fm-j + 2\ 
- 2 ^ (m-j + 2)2— V m-j + 1 A J / 

Thus formula f lTOj) gives explicit solution of equations (I2])-(|4]) for the distribution rijS for any 
value of 9. 

Using known assymptotics for Motzkin numbers (after Benoit Cloitre, see jsl) 

c(z) ~ v/374A^ (23) 

we can present explicit formula for asymptotic behaviour for the rii distribution. In the limit 
case 9 — > 0, while N9 = const, 

= jc^-i {N9){^yc,.u (24) 

and using ( 123|1 one obtains 

rii+i ~ ^. (25) 

Thus we found one more example of inverse-power distribution with the power equal to |. 
Author would like to thank prof. Zbigniew Czechowski for fruitful discussions. 
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